Modified Kuramoto-Sivashinsky equation: stability of stationary solutions and the 

consequent dynamics 
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We study the effect of a higher-order nonlinearity in the standard Kuramoto-Sivashinsky equation: 
d x G(H x ). We find that the stability of steady states depends on dv/dq, the derivative of the interface 
velocity on the wavevector q of the steady state. If the standard nonlinearity vanishes, coarsening 
is possible, in principle, only if G is an odd function of H x . In this case, the equation falls in the 
category of the generalized Cahn-Hilliard equation, whose dynamical behavior was recently studied 
by the same authors. Instead, if G is an even function of H x , we show that steady-state solutions 
are not permissible. 

PACS numbers: 05.45.-a, 82.40.Ck, 02.30.Jr 
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I. INTRODUCTION 

One of the most prominent and generic equations 
that arises in nonequilibrium systems is the Kuramoto- 
Sivashinsky (KS) 0,110 



H t + H x 



equation: 
H xx + HH X = 0, 



(1) 



where H is some scalar function (like the slope of a one- 
dimensional growing front), and differentiations are sub- 
scripted. The linear stability analysis of the KS equa- 
tion (by looking for solutions in the form of e l i x + ult ^ 
yields u = q 2 — q 4 . The (linearly) fastest growing mode 
has a wavenumber given by q u = l/\/2 (obtained from 
dgUJ = 0). For a large box size the KS equation is known 
to exhibit spatiotcmporal chaos. The chaotic pattern sta- 
tistically selects a length scale which is close to 2Ti/q u : 
in fact, the structure factor (\H q \ 2 ), where (•••) desig- 
nates the average over many runs and H q is the Fourier 
transform of H, exhibits a maximum around q = q u . 
Other nonequlibrium equations are known, however, to 
exhibit different dynamical behaviors: just to limit to 
one-dimensional systems, we may have coarsening, a di- 
verging amplitude with a fixed wavelength, a frozen pat- 
tern, travelling waves, and so on. 

An important issue is the recognition of general crite- 
ria that enable to predict whether or not coarsening takes 
place within a class of nonlinear equations, without hav- 
ing to resort to a forward time dependent calculation. In 
recent works 0, 0, @ we have considered several classes 
of one-dimensional Partial Differential Equations (PDE), 
having the form H t = M[H], where M is a nonlinear op- 
erator acting on the spatial variable x. 

Sometimes, even in the presence of strong nonlinear- 
ities, the search for steady states reduces to solving a 
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Newton-type equation, H xx + V(H) = 0, where H is 
some function of H. In these cases, X(A), giving the 
dependence of the wavelength A of the steady state on 
its amplitude A, is a one- value function (see Fig. 2] full 
lines), and the criterion for the existence of coarsening is 
expressed in terms of the derivative X'(A). 

It has been shown that A' (A) has minus the sign of 
the phase diffusion equation. Thus X'(A) > corre- 
sponds to a branch which is unstable with respect to 
the phase of the pattern, entailing thus coarsening. The 
situation is more complicated when X(A) exhibits a fold 
(see Fig. n dashed line). This event occurs, e.g., in the 
KS equation and in the Swift-Hohenberg equation. As 
for the KS equation, which is the topic of this paper, 
Nepomnyashchii pj has shown that the forearm part of 
the curve X(A) with positive slope, A' (A) > 0, is an un- 
stable branch. This result holds only for the pure KS 
equation, however. 

The aim of this paper is the following, (i) Firstly we 
shall extend the result of Nepomnyashchii 1] to a gen- 
eralized form of the KS equation, which includes higher 
order nonlincarities. As a way of example, the next lead- 
ingterm in the KS equation (H X H XX ) has been analyzed 
in [2j, and it has been shown that this term significantly 
affect dynamics; for example the profile may exhibit deep 
grooves. We shall consider a more general form of the 
modified KS equation by adding a term like d x G(H x ) 
(this includes as a particular case the term H X H XX ). We 
find that the stability of the steady state solutions de- 
pends on v'(q), where v is the average interface velocity, 
(ii) If the standard KS nonlinearity vanishes (no HH X 
term), then there is coarsening if G is an odd function. 
In this case a mapping of the equation onto a generalized 
Cahn-Hilliard equation is straightforward, (iii) If G is 
even (still in the absence of the standard nonlinearity), 
we show that there exists no steady-state periodic solu- 
tion, as attested by numerical simulations for G = H x . 
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FIG. 1: The periodic stationary solutions of some classes of 
nonlinear equations satisfy a Newton-type equation, H xx + 
V(H) = 0, and the resulting X(A) is a single-value function. 
According to the explicit form of V(H), we may have the 
different curves shown as full lines. The dashed line, which 
displays a fold, corresponds to a non single- value function. 
It comes out in the Kuramoto-Sivashinsky and in the Swift- 
Hohenberg equations. Units on both axes are arbitrary. 



II. THE MODIFIED KS EQUATION 
A. The method 

We study the following equation: 

H t + c 4 H xxxx + c 2 H xx + aHH x + (3d x G(H x ) = 0, (2) 

which reduces to the standard Kuramoto-Sivashinsky 
equation when G = 0. A rescaling of x, t and H always 
allows one to reduce the equation to a one parameter 
equation, which can be absorbed into a redifmition of G. 
However, for the sake of clarity, we do not get rid of a, 
so that we write 



Ht + H x 



H xx + aHH x + d x G(H x ) = 0. (3) 



It is also useful to rewrite Eq. J3J using the variable u, 
with u T = H: 



u t + u xxxx + u xx + %u x + G{u xx ) = 



(4) 



where the integration constant vq can be canceled out by 
the transformation u(x, t) — > u(x, t) — vot. 

Within the H— formulation, the average velocity 
d(H)/dt vanishes, because Eq. @ has the conserved 
form H t = —d x {. . . ). Within the it— formulation, 



d(u) 
dt 



(ul) + (G(u xx )) 



(5) 



We start from a stationary solution of period q, H(x), 
and perturb it with adding h(x) exp(— cut). The function 
h{x) therefore satisfies the linear equation 

+ h xx + ot{hH) x + (h x G(H x )). x = uih, (6) 



where G = G' and whose coefficients are periodic with 
period A = 27r/g. According to the Floquet-Bloch the- 
orem, the solution has the form h(x) = exp(iKx)F(x), 
where F(x) has the same period as H{x). 

We are interested in weak modulations of long wave- 
length, i.e. with K <C q. It is therefore convenient to in- 
troduce the reduced wavevector Q = K/q and the phase 
4> = qx. The equation determining the steady state H(<j>) 
is 

<J 4 i^00 + q 2 H 4>4> + aqHHtf, + q 2 H 4>4> G(qH, p ) = (7) 

and Eq. © reads 

q A h (t>(j)< j >( j, + q 2 h cj}t j ) + qa^hH)^ + q 2 (h (j> G(qH^)),j ) 

= Ch = ujh, (8) 

with h = exp(iQcj))F((j)). 

The final step is to expand both lj and F((f>) in powers 
of Q, 

F(<t>) = F Q {0)+QF 1 {0)+Q 2 F 2 ( ( j ) ) + ... (9) 
CJ = UJq + Qwi + Q 2 oj 2 + ■■■ (10) 

and to solve Eq. JSJl at the lowest orders in Q. 



B. Zero order 

The differential equation determining the steady state 
H{x) does not depend explicitely on x, so that H(x + xq) 
is a solution as well. This symmetry implies that Ch = 
ujh (see Eq. (JSJ) is solved by h = and co = 0, as can 
be easily checked by taking the <j> derivative of Eq. J7J). 
Since the zero order equation is simply 



CF Q = 0, (11) 



we obtain Fn — Ha 



C. First order 

The equation for F\ reads 

CFt = (12) 
UlH' - iq[4q 3 H"" + 2q(l + G)H" + {aH + q 2 H"G')H'], 

where we have used the shorthands H 1 = , H" = 
H^, ... . If we differentiate Eq. Q with respect to q, 
we get a similar equation, 

CH q = -(4:q 3 H"" + 2qH" + aHH' + 3(3q 2 H'H"). (13) 



The comparison of Eqs. (|12I13|1 suggests to look for F± 
under the form F\ = iqH q + c, where c is a constant. We 
easily find c = u>i/(aq), so that 



Fx 



aq 



iqH q 



(14) 
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This result shows that in the absence of the standard 
nonlinearity, a — 0, u>i should vanish whatever G is. 

We might have started with an even more general 
Eq. 0, replacing the standard nonlinearity (aHH x ) = 
d x (§H 2 ) with d x P{H). The term aHH' in the right 
hand sides of Eqs. I|12|) and i|13|) would be replaced by 
P(H)H' , with P = P' . However, in the general case P' 
is not a constant: therefore, it is not possible to look for 
a solution F% — iqH q + c. 

D. Second order 



where 

Cu = u xxxx + u xx + a(Hu) x + [G{H')u x ] x (20) 

and 

£t u = u xxxx + u xx - aHu x + [G{H')u x ] x . (21) 

The equation Cu = is solved by u = H 1 , but we do 
not know the solution of C^u = 0. & ^ C because of the 
a-term. If a = 0, such term is absent and C) = L. This 
case is treated in the next subsection. 



The equation for F 2 has the form 

CF 2 +q 2 {l + G)(-F + 2iF[) + (aqH + g 3 #"G')iiq 
= u 2 F +u 1 F l , (15) 

r 

LF 2 = lu 2 H' + iquj 1 H q + q 2 H' + 2q 3 H' q -iuj 1 H 
-i{uo 1 /a)q 2 H"G'{qH') 
Wj{aq) + q 2 G(qH')H' + 2q 3 G(qH')H' 



+aq 2 HH q + q 4 H"G'(qH')H q . 



(16) 



Now we take the 2-7r-average of the previous equation, 
getting 



-± + q 2 (G(qH')H') + 2q 3 (G(q(H')H' q ) 
+aq 2 (HH q ) + q A {H"G'{qH')H q ) = 0. 



(17) 



Finally, we obtain 
d 



-aq 3 — 
aq 



2 X 1 



(G(qH')) =aq 3 -v. (18) 
aq 



This result proves that u>x = if a = 0, whatever the 
function G is. Consider the case a/0 (we are at liberty 
to choose a > 0). Since uj < signals an instability, one 
sees that if j^v > 0, the periodic solution is unstable, 
because there is a real solution lo\ < 0. This generalizes 
the result of obtained for the pure KS equation, to the 
higher order KS equation. It is only in the pure KS limit 
that the spectrum of stability is related to the slope of the 
steady amplitude. In the higher order equation, however, 
this ceases to be the case. Instead we should replace the 
amplitude by the drift velocity, a quantity which can be 
still obtained from pure steady-state considerations. 



E. Determination of u>2 

The determination of uj 2 implies the resolution of the 
differential equation 



C)u = 0, 



F. The case a = 

If a = 0, stationary solutions are determined by the 
equation 



G{H X ) = C, 
where C is a constant. If h = H x , then 
h xx = —h — G{h) + Cx. 



(22) 



(23) 



The equation for H (x) admits periodic solutions only 
if we take C = and if G(h) is an odd function, so that 
h(x) itself is periodic and has zero average for any initial 
condition (such that the solution is bounded). 

In the same limit a — 0, the full PDE JSJ) writes 



H t = —d x [H xxx + H x + G(H X )] 



(24) 



and taking the spatial derivative of both terms, we get 



h t = -d xx [h xx + h + G(h)}, 



(25) 



where h = H x . We have therefore got a generalized 
Cahn-Hilliard equation, whose dynamical behavior is 
known to show coarsening if and only if the wavelength A 
of steady states is an increasing function of their ampli- 
tude A |5j . We have reobtained the same result following 
the method discussed in this Section. 



III. FINAL REMARKS 



(19) 



Present and recent work 4, 5, 6] has the main objective 
to find general criteria to understand and anticipate the 
dynamics of nonlinear systems by the analysis of steady 
state solutions only. For some important classes of PDE, 
which have the common feature of a single-value X(A) 
function, the criterion is based on the sign of the deriva- 
tive A'(j4). In this short note we have considered a mod- 
ified, generalized Kuramoto-Sivashisnky equation, where 
the A (A) curve is not single- value, because it displays a 
fold. We have therefore established a different criterion, 
based on dv/dq, the derivative of the interface velocity 
on the wavevector q of the steady state. 
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If the standard KS nonlinearity is absent we can say 
more. If the nonlinearity d x G(H x ) corresponds to an 
even function G, the equation does not support peri- 
odic stationary solutions and this prevents coarsening in 
principle: we have a pattern of fixed wavelength and di- 
verging amplitude. Instead, if G is an odd function, the 



equation falls in a previous studied class, the generalized 
Cahn-Hilliard equation, which can show different behav- 
iors according to the form of X(A). Finally, it is an im- 
portant task for future investigations to see whether or 
not information of the types presented here and in 0, Q 
have analogues in higher dimensions. 
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